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Abstract 

We measure the displacement response of a frictional multicontact inter- 
face between identical polymer glasses to a biased shear force oscillation. We 
evidence the existence, for maximum forces close below the nominal static 
threshold, of a jamming creep regime governed by an ageing-rejuvenation 
competition acting within the micrometer-sized contacting asperities. Quan- 
titative analysis of the creep curves suggests that another such mechanism 

might be at work within the nanometer-thick adhesive junctions. 
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Solid friction between two macroscopic solids is commonly characterized 
in terms of (i) a static force threshold below which no relative displacement 
is supposed to take place, (ii) a dynamic friction coefficient, measured in sta- 
tionary motion. However, recent experiments performed on multicontact in- 
terfaces (MCI), i.e. interfaces between two solids with rough surfaces, pressed 
together under a load N, have revealed that ||: 

(i) for shear forces F such that F/N <C fi s , where pL s is the static threshold, 
the pinned interface responds elastically, via the reversible deformation of the 
contacting asperities, 

(ii) for F < fJ> s N, creeplike irreversible sliding is observed. 

Clearly, the study of the latter regime of incipient sliding should give access 
to precise information about the underlying pinning/depinning dynamics. 

It is now well established § that, for a MCI, the variations of the friction 
force F = o" s S r (with E r the real area of contact, and a s the interfacial shear 
stress) are governed, at low velocities, by the competition between two effects: 
(i) an age strengthening effect resulting from the logarithmic creep growth, 
under the high (geometry-enhanced) normal load, of the sparse contacts be- 
tween load bearing asperities. When motion starts, contacts get gradually 
destroyed, after a lifetime or age <3?, and replaced by fresh ones. So, while 
the interface sits still, it ages (strengthens), when it slides, it rejuvenates 
(weakens). Full refreshment occurs, on average, after sliding a micromet- 
ric memory length Dq. (ii) Weakening when sliding is counteracted by the 
velocity-strengthening interface rheology: a s {x) = a s o [1 + a\n(x/Vo)], with 
x the instantaneous sliding speed. It results from thermally activated prema- 
ture depinning events within the nanometer-thick adhesive junctions between 
contacting asperities. 

Both effects yield logarithmic variations of F. One thus expects creep to 
exhibit a strong, exponential, sensitivity to forces close to the nominal static 
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threshold. 

Now, the experiments reported in Q were performed under static loading 
through a spring of finite stiffness. As such, they did not provide a force 
control fine enough to study incipient creep accurately. So, we have chosen, in 
this work, to probe it via the response to a biased a.c. shear force. With a bias 
Fdc ^ ftsN, and an amplitude such that the maximum force F max lies in the 
tangential creep range, the interface should experience, during each oscillation 
period, an alternation of two regimes: (i) for F close to F max , a sliding 
phase during which rejuvenation is at work, yielding an age variation A^ s ud e , 
and (ii) as F decreases, the slip velocity decreases quasi-exponentially, and 
the slider enters, for the rest of the period, a quasi static phase where age 
grows linearly with time by an amount A$> s t a t- Such a competition between 
rejuvenation and ageing is akin to that invoked to model soft glassy rheology 
(SGR) §. 

If, say, s nde + A<3? stat < 0, the interface will weaken, leading to a larger 
slip during the next period, etc. One thus expects the dynamics to bifur- 
cate between self- accelerated unlimited slip and self-decelerated, jamming ||, 
creep. The experiments reported below fully confirm this qualitative scenario. 
Moreover, quantitative analysis of the data, based upon the Rice-Ruina (RR) 
model ||, allows us to show that the rejuvenation-ageing process cannot be 
fully ascribed to variations of S r , leading us to conclude that the interfacial 
rheology is most likely, itself, of the SGR type. 

Experiments. — The experimental setup has been fully described in ||; it 
is sketched in the inset of Figjlja. Two poly(methyl methacrylate) (PMMA) 
samples, with lapped surfaces of roughness R a = 1 /xm, are glued on a slider 
and a track and form the multicontact interface. The slider, of nominal area 
4 cm 2 , rests on the track, inclined at 9 = 20° from the horizontal. The 
tangential (Fdc) to normal (N) load ratio "fdc = tan 9 = 0.36 is well below the 
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static threshold fi s ~ 0.6 (see [||) and no sliding occurs. Imposing a harmonic 
motion to the track then results in an inertial shear loading of the slider, of 
amplitude F ac = ^f ac N, with j ac < 0.5. The frequency / = 80 Hz is chosen 
well below the natural frequency of the slider-interface system, /o = 800 Hz, 
so that the inertia associated to its relative motion in the track frame can 
be neglected. We measure the displacement X of the center of mass of the 
slider by means of a capacitive displacement gauge, with a noise amplitude 1 
nm over its whole 0-500 Hz bandwidth. In order to prepare the system in as 
reproducible as possible an initial state, the slider is placed on the track and 
a large j ac is then imposed, in order to make it slide a few micrometers in the 
direction of Ffa. The harmonic force is then suddenly stopped, which results 
in an elastic recoil of the contacting asperities 0. This method reduces the 
relative dispersion on interfacial stiffness values (see ||) to only 10%. This 
value agrees with the expected statistical dispersion due to the finite number 
of load bearing contacts, which we can estimate to be of order 50. A time 
twait is then waited before reswitching the harmonic shear loading, either as 
a linear ramp of amplitude, until gross sliding occurs, or as a step with rising 
time < 0.1 s. 

Results. — The displacement response X(t) of the slider to a ramp 7 ac (i) is 
illustrated on Fig.||a. Also shown is the average displacement X^ c measured 
by filtering X(t) through a low-pass filter of cutoff frequency 8 Hz. In region 
(I), the slider oscillates about a constant average position, no irreversible slip 
occurs: the MCI responds elastically. Region (III) corresponds to accelerated 
sliding. In region (II) between these two regimes, the average displacement 
Xdc, i-e. the slipped distance, increases continuously. We use this ramp test 
to define a threshold 7 S = (i ? d c + F ac )/N such that the average sliding velocity 
dXd c /dt = 100/™i.s _1 . We thus obtain, for t wa u = 600 s and 7 ac = 0.1 s , 
7 S = 0.59 ± 0.03. The scattering, of order 10%, is consistent with that of the 
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stiffness. 

One can see, on Fig.[I]a, that the intermediate creep regime (II) corresponds 
to a narrow range of 7 max = 7<j c + 7 ac . Slow creep is studied by choosing a 
value ^fmax in this range, setting, at t = t wa u, the amplitude stepwise to 7 ac , 
and recording X^cit). 

The set of creep curves displayed on Fig.||b all correspond to t wa u = 300 s, 
7max = 0.54. The large dispersion between X<± c for various runs must therefore 
result from the statistical dispersion of the MCI initial state. After slipping 
by a finite amount, of order 10-100 nm, over the rising time (0.1 s) of j ac , 
the slider performs a slowly self decelerating creep. After, typically, 10 4 s, the 
slip velocity has decreased to non-measurable values, indicating a saturating, 
jamming dynamics. We attribute the large dispersion of the creep curves 
to the expected above-mentioned exponential sensitivity of the dynamics to 
7s — Imax- A direct confirmation of this is obtained from the experiment 
presented on Figjljc: it shows that a 3% step of ^ ma x turns quasi-jamming 
into accelerated sliding. 

These ideas can be checked in a more quantitative way as follows: as long 
as the age of the MCI has not been appreciably modified by the creeping 
dynamics itself, we expect the characteristic time for creep, t c , to be that 
for thermally activated depinning of a typical nm 3 pinned unit within the 
adhesive layer: ln(t c ) ~ Cst + AE/kT, where AE is the energy barrier to 
be jumped by an element under reduced load "fmax- Close to the depinning 
threshold § AE/kT ~ ( 7s 

— ^ max) I A. The RR rate parameter A — acr S Q/p, 
with p the average pressure on the microcontacts, has been measured, for 
PMMA/PMMA, to be 0.013 |§. Hence, we expect the slipped distance in run 
(i) to scale as X l dc (t) = X(t/d), where ln(Ci) ~ (7* - jmax)/A. Fig.| shows 
the set of creep curves resulting from such a scaling. Indeed, the collapse onto 
a master curve is very good in the short time range. Moreover, we find the 
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maximum spread of the scaling factors max | m(Cj/C,-)| ~ 4.4, in excellent 
agreement with max|A(7* — ^max)\/A ~ 4.6. 

Discussion. — We now analyze our results within the framework of the 
Rice-Ruina phenomenology ||, which has proved to account very well for 
the stick-slip frictional dynamics of MCIs Q. It models the above described 
rejuvenation- ageing process and the velocity strengthening interface rheology 
as follows. 

The friction coefficient reads: 

fj, = F/N = fi Q + A In (£) + J Bln^ <I>i " 



V J VA) 

is the interface age, no the friction coefficient at reference velocity Vo- The 
instantaneous interfacial sliding velocity x is related to the center of mass 
position by x = d(X — F/K)/dt, with k the interfacial elastic stiffness Q. 
The age <& evolves according to: 

^ -, 

<3? = 1 

D 

On the r.h.s. of Eq.(|2|), the first and second terms correspond, respectively, 
to time ageing and slip rejuvenation. 

We have performed numerical integrations of this set of differential equa- 
tions to calculate the slipped distance X^ c (t) with, in Eq.([T]), F/N = 
Idc + 7acCos(o;t), and initial conditions for slip and age Xd c (0) = 0, $(0) = 
twait = 300 s. We have used for the memory length Dq and the dimensionless 
parameters A and B the values Dq = 0.42 ^m, A = 0.013, B = 0.026, obtained 
from previous measurements on PMMA |J. We choose Vq = l/im.s" 1 . Due 
to the exponential amplification by the creep dynamics of the small variations 
of the absolute friction level between various runs, must be left free. This 
unique fitting parameter is tuned so as to adjust the calculated and measured 
values of X^c at the end of the run. A typical example of such fits is shown 
on Fig.||. It yields //o = 0.42865 |pX|| , fully compatible with previous data [0. 
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These fits appear to be very good at long times t > 1000 s, that is in 
the quasi-jammed regime where static ageing becomes dominant. The corre- 
sponding asymptotic dynamics can be analyzed directly. From Eq.(|l|): 



x 



D ^ 
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7 ac (cos(q;£) - 1) 
A 



with [3 = B/A. 

Let us consider the oscillation period centered at t±, with j(t%) = ^ m ax- 
The velocity x only takes on significant values for 7 ~ Jmaxi *.e. f° r l T l = 
\t - ti| < T. So, in Eq.@, cos (art) - 1 « -u 2 t 2 /2, and: 
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with the constant $f = (Do/Vq)^ exp [(j ma x - Ho) /A] and wr c = (2A/7 ac ) 1 / 2 . 

To lowest order, the slip-induced variation can be neglected, and Eq.(|2|) 
yields: 



where the integration constant is the age at some "initial" time within the 
quasi-jammed regime. 

From Eqs.(||) and (||), the increment of slip over this period: 

J2- 
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where the integration can be extended to infinity due to the gaussian decay 
of x. Then, for t\ S> 0: 
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From this, the sliding velocity Xd c , coarse-grained over the period T: 



and we obtain for the slipped distance: 



Since, for our system, (3 = 2, we thus expect the creeped distance to approach 
its saturation level as 1/t. 

Our experimental data are seen on the insert of Fig.|| to fit the jamming 
asymptotics predicted by the RR model with excellent accuracy. 

However, for all the experimental runs, we find (see Figj3|) that, although 
the overall shape of Xdc(t) is reasonably well described by the RR fits, the 
agreement is clearly not quantitative at short times (t < 1000 s). In this 
time bracket, the RR model systematically underestimates Xd c (t), hence the 
rejuvenation efficiency of slip. In particular, it by no means accounts for 
the fast increase of X^c, which amounts typically to ~ 10% of the total slip, 
occurring over the stepping time. 

This strongly hints at the fact that the RR model, while it very well de- 
scribes established sliding, misses some important feature of incipient sliding. 
We suspect that this missing feature might be slip-induced rejuvenation, i.e. 
dynamical weakening of a s , within the nanometer-thick adhesive junctions 
themselves. Indeed, these certainly have an amorphous solid structure when 
pinned, and they flow beyond a stress threshold. As such, they can reason- 
ably be expected to behave as soft glassy materials, whose rheology is now 
interpreted || in terms of structural ageing/rejuvenation competition. If this 
turns out to be the case, two such mechanisms would be at work in the MCI 
solid friction, on the two scales of, respectively, the micrometric asperities, 
and the nanometric pinning units. 

This issue is, in particular, of primary relevance to the modelling of the 



dynamics of interfacial shear fracture [11,12] . We plan to investigate it by 
studying the frictional dynamics of rough PMMA sliding over smooth hard 
glass. Since, with such a system, the microcontact population is unaffected 
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by motion, rejuvenation, if observed, will have to originate from the adhesive 
junctions. 
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FIGURES 

FIG. 1. (a) Instantaneous (line) and averaged (•) displacement response of the slider to a biased 
oscillating shear force of ramped reduced amplitude 7ac(i) (°)- The bias "fdc = 0.36. The arrow 
indicates the point at which the averaged velocity reaches 100 ^m.s -1 . Inset: experimental setup, 
(b) Recordings of the average creep displacement X^ c for four runs performed under identical 
nominal conditions: 7^ c = 0.36, "y ac = 0.18, t wa u = 300 s. The wide scattering of the curves results 
from dynamical amplification of the statistical dispersion of the interfacial strength, (c) Transition 
from jamming to unbounded slip (•) triggered by a 3% jump of ^ ma x = 7tfc + lac (°)- 

FIG. 2. Scaled plot of 7 creep curves (same nominal conditions as for Fig. lb). The reference 
run (C = 1) corresponds to (o) symbols. 

FIG. 3. An experimental creep curve (•) and its fit according to RR model (line). With 
7d c = 0.36, 7 ac = 0.18, t wa it = 300 s. Inset: the same plotted versus 1/t, the dashed line indicates 
the quasi-jamming asymptotics (see text). 
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